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Gas lenses are being considered as focusing elements of beam-waveguides. 
Since very many lenses are needed to form a long waveguide, it is reasonable 
to consider periodic arrangements of gas lenses. Such periodic structures 
might operate with a gas stream which flows through all of the lenses in 
succession. A periodic temperature distribution in the gas results which is 
different from that of single gas lenses considered in two earlier papers. 3 A 

This paper analyses the ray optics properties, such as focal length and 
principal surfaces of the gas lenses, of two types of alternating gradient 
focusing systems. One system consists simply of a succession of hot and 
cold tubes. The other system results from the first by insertion of heat in- 
sulating tubes of equal length between the hot and cold tubes. 

I. INTRODUCTION 

The beam-waveguide described by Goubau 1 appears as a promising 
device to transmit light over long distances. However, to reduce the 
power loss due to absorption and reflection, which is inevitable with 
lenses made of solid dielectrics, gas lenses have been proposed 2 3 instead 
of the solid lenses used by Goubau. 

Two earlier papers 3,4 discussed the properties of a particular type of 
gas lens. This tubular gas lens consists of a warm tube into which a 
cooler gas is blown. The thermal gradients in the gas lead to density 
gradients which give the structure the properties of a positive lens. 

Ref . 4 discusses the focal length and principal surface of this gas lens 
for the case that the gas enters the lens at a constant temperature. It 
was shown that this device, when operated under optimum conditions, 
acted as an optically rather thin lens with moderate lens distortions. 

The present paper extends the earlier analysis in several ways. We 
consider periodic lens structures. Such a structure results if hot and cold 
tubes are alternated to form a long, periodic structure. The gas is heated 
and cooled periodically giving rise to periodically arranged positive and 
negative lenses. A periodic structure of this type represents an alter- 
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nating gradient focusing system. 5 In general, the temperature of the gas 
entering the hot or cold tubes will not be constant over the cross section 
of the tube so that the earlier results are no longer applicable. The gas 
temperature in the periodic structure will also be periodic and will de- 
pend on the temperatures of the hot and cold tubes as well as on the flow 
velocity. 

It is our aim to compute the temperature distribution in such periodic 
structures and use it to determine the properties of the equivalent 
lenses which describe the ray optics of the alternating gradient focusing 
systems. 

The equivalent lenses are rather complicated. They are neither op- 
tically thin nor free of distortions. Further investigations are required 
to determine the guidance properties of an alternating gradient focusing 
system with imperfect lenses of this type. 

We discuss two types of periodic gas lens systems. In one case we 
assume that hot and cold tubes of equal length are directly adjacent to 
each other. The other type is an alternating gradient beam-waveguide 
which consists of hot and cold tubes which are separated by tube sec- 
tions made of an ideally heat insulating material. For simplicity it is 
assumed that the insulating sections are as long as the hot and cold 
tubes. The assumption of a perfectly heat insulating material is an over- 
idealization since hardly any material conducts heat more poorly than 
gases. It is intended as an approximation to the real situation of im- 
perfect heat insulators. 

In the insulating tube sections, the gas has a chance to equalize its 
temperature. As it does so rather rapidly, we again have the case of hot 
and cold tubes being fed by an input gas at a constant temperature. 
However, the insulating sections act also as lenses in the same sense as 
the hot or cold tubes by which they are preceeded. Therefore, it is not 
surprising that some improvement of efficiency results if heat insulating 
tube sections are used to separate the hot and cold tubes. But this ad- 
vantage is not very striking; and, since this analysis assumes ideally 
insulating tubes, it is not certain how much of a real advantage can be 
gained by using this construction. Considerably more experience is 
needed before a decision can be made. 

This analysis again neglects all convection effects in the gas lenses. 

To be able to distinguish which of the two structures is being dis- 
cussed we will call the structure using hot and cold tubes without heat 
insulating tube sections the simple periodic structure while the second 
case which includes insulating sections will be called the extended periodic 
structure. 
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II. RAY TRACING 



Before entering into a discussion of the simple and extended periodic 
structures, the ray tracing technique used to determine the focal length 
and principal surface of the lenses will be explained. 

The trajectory of the light ray in the gas lens is given by the ray 
equation 6 



- (n -\ 

(Is \ dsj 



— J = grad n ( 1 ) 



r = the position vector leading from an arbitrary origin to points 

on the ray. 
n = index of refraction. 
s = length coordinate measured along the ray. 

We limit ourselves to rays which are very nearly parallel to the axis 
of the structure which is used as the ^-coordinate so that we can replace 
s by z.* Assuming angular symmetry it is sufficient to consider the vector 
component in radial direction /• perpendicular to the z-axis. Finally, we 
neglect the term 

(dn/dz) (dr/dz) 

because dr/dz « 1 for rays which are nearly parallel to the z-axis and 
also because the variation of n in the z-direction will generally be smaller 
than that in the r-direction. 

(dn/dz) « (dn/dr). 

With these assumptions we obtain from (1) 

^' = i^. (2) 

dz- n or 

However, since we are only interested in gases where n — 1 « 1 we can 
safely write 

d*r _ dn /,,, 

dz 2 dr 

The index of refraction depends on temperature in the following way: 

n - 1 = (n -Dp (4) 

To is the absolute temperature at which n is measured while T is the 

* The error caused by this approximation is estimated in Ref. 4. 
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absolute temperature for which we want to determine n. It follows 
from (4) that 

dn i ,\ To dT 

_- -(n a -l) ¥2 -. 

The value of the absolute temperature varies only slightly through- 
out the gas so that T can be replaced by a suitable average temperature. 
It is convenient to choose T equal to this average temperature which 
should be chosen as 

T Q = $(T h +T„) (5) 

Th = temperature of hot tubes 
T c = temperature of cold tubes 

This leads us to the ray equation 

ttr = wo- 1 dT (6) 

dz* T dr ' 

Equation (6) is our starting point for the ray optics of the gas lenses. 
Since dT/dr is a complicated function of r and z, it is difficult to solve 
(6) analytically so that we content ourselves with numerical solutions 
obtained by means of an electronic computer. 

Rather than expressing our results as functions of z, we want to ob- 
tain them as functions of the on-axis gas velocity v normalized by a 
suitable constant V. (This representation was also used in Refs. 3 and 
4.) We define V (L) by 

v /V(L) = a/oL (7) 



with 



a = k/avopCp (8) 



a 



= tube radius 
L = tube length 
k = heat conductivity of the gas 
p = (average) gas density 
c p = specific heat at constant gas pressure. 

Equation (7) shows that v /V(L) is inversely proportional to the length 
of the tube. Therefore, it is convenient to introduce a variable 

u(z) = alaz (9) 



which at z = L equals 

Using 

equation (6) becomes 
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w-.a)-jfc. do: 



x = r/a (11) 



d 2 x 2 dx _ _ 1 n - 1 dT 

dv? + udu u 4 a 2 T dx' K ' 

Equation (12) is used to obtain x and dxldu as a function of W from 
which the focal length / and principal surface p can be computed. 

As shown in Fig. 1 (a), we follow the ray from a point z = Z\ to z = 
«2 , corresponding to u = U\ and u = Ui , through the tube anticipating 
the case of lenses which are not bounded by planes through the ends of 
the tube but by surfaces inside of the tube to be defined later. 

The definition of focal length and principal surface can be seen from 
Fig. 1 (a). The principal surface is obtained by following a ray, which is 
incident parallel to the axis (dx/dz = at 2 = Zi), through the lens. 
If we extend the direction of the ray entering the lens and the direction 
of the ray leaving the lens at z = z 2 by straight lines back into the lens 
we obtain a cross-over point which defines a point on the principal 
surface. The distance p of this point measured from the beginning of the 
tube as a function of xi , the input position of the ray, describes the 
principal surface. The distance p+ for rays traveling in the same direc- 
tion as the gas flow is given by (Fig. 1 (a)) 

x(z 2 ) - s(z i) 

p+ = 2 ' +L ~^A (13) 

dz),=,, 

or, expressed in terms of it and W rather than L 

ft .»■ + 1 + '<»>-'<*> y. 

L L ■2{dx\ (14) 



du/ u =u 3 

We define the focal length as the distance from the intersection of the 
incoming ray (extended in a straight line) with the principal surface 
to the point at which the outgoing ray crosses the axis of the structure. 
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(b) 

Fig. 1 — Geometry of a single gas lens showing the definition of focal length 
and principal surfaces. 



Therefore, we have 



/+ = - 



aQi) 
(dx\ 
\dz)z = . 



(15) 



or, in terms of u and W, 



L 2 



x(ui) 



U 2 \ -=- 



(16) 



(dx\ ■ 

\du/u=u 2 

Similarly, we obtain from Fig. 1(b) the principal surface and focal 
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length for the ray traveling in opposite direction to the gas flow 

t - 5 - * 



x(u 2 ) - x(ui) 



/dx\ 

\du/ u=lll 



L L u s (^ (17) 



and 



f- = _ w *M 

L ' u 2 (dx\ IS 






The principal surfaces p + and p_ do not, in general, coincide. If they 
are identical the lens is called optically thin. The separation between 
the two principal surfaces is an indication of the optical thickness of the 
lens. 

III. THE SIMPLE PERIODIC STRUCTURE 

3.1 Temperature Distribution 

The simple periodic structure consists of alternating hot and cold 
tubes (Fig. 2). In order to compute the equivalent positive and negative 
lenses of this structure we first have to determine the temperature 
distribution. 

The temperature distribution is given by a series expansion 7 which, 
in the hot tube, reads 

00 

Tdx,u) = T h - £ AnR»(x) exp (-0„7w) (19) 

X > U > W = =5- V ° 



V(L) 
with T h being the wall temperature of the hot tube and u = a/oz; and 



GAS 
FLOW 



Z=0 Z = L Z = 2L Z = 3L Z = 4L 

Fig. 2 — Sequence of hot and cold tubes comprising the "simple periodic 
structure". 
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in the cold tube 

T*(x,u) = T e - £ BAW exp |- jB. f Yi - iM (20) 

with 

IF > u > £W 

T c = wall temperature of cold tube 
W = a/oL. 

The functions R n and the eigenvalues /3„ are discussed in the appendix. 
The coefficients A„ and B„ have to be determined so that the temperature 
is a periodic function in the simple periodic structure of Fig. 2. 

To simplify the determination of A„ and B n we limit ourselves to 
sufficiently long tubes or slow enough flow velocities so that the first 
term of the series expansions (19) and (20) are sufficient to describe 
the temperature distribution at the end of each tube accurately. This 
condition is expressed by the requirements 

w = vh <10 - (21) 

The periodicity condition requires that the temperature at the end of 
the cold tube (z = 2L or u = \W) equals the temperature at the be- 
ginning of the hot tube (z = or m = =» ) 

T x {x,«>) = TtMW)- . (22) 

In addition, we have to require that the gas temperature passes con- 
tinuously from the hot to the cold tube 

r,(.r,TF) = T 2 (x,W). (23) 

The conditions (21) to (23) allow the determination of the constants 
A n and B n 

2(T h - T e ) f, 5o 

30 A=i 
and 

B n = - A n (25) 

with 

1 n = 




f*R*\ ■ —(** ' V '' 



■ " 10 n * 0. 
The eigenvalues /3„ and the derivation dR /d(3 are given in the appendix. 
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The temperature distribution in the hot tube is shown in Fig. 3(a) 
for v /V(L) = 5 and in Fig. 3(b) for v /V(L) = 10 as a function of 
normalized radium x = r/a. The various curves in each figure correspond 
to different positions along the tube axis. At z = the temperature 
distribution is identical to that at the end of the cold tube. The tem- 
perature is cold on the wall and warmer in the center of the tube. As 
we follow the temperature distribution deeper into the hot tube we see 
that the temperature on the wall changes instantly from its value equal 
to the wall temperature of the cold tube to that of the hot tube. How- 
ever, the slope of the temperature distribution close to the tube axis 
remains negative tor quite some length. This means that the gas in the 
hot tube acts like a negative lens close to the input end of the tube. It 
takes some distance to reverse the negative temperature gradient which 
the cold tube imparted to the gas. In fact, there exists a neutral surface 
in the hot as well as the cold tube which is denned by the points where 
the temperature gradient dT/d.r = 0. On this surface the gas acts neither 
as a positive nor negative lens. The neutral surface separates the region 



X 

Th 




Fig. 3(h) -Temperature distribution in the simple periodic structure as a 
function of normalized radius x = r/a at various cross sections z/L in the tube. 
The normalized flow velocity is v /V(L) = 5 and (T h - T r )/T k = 0.155. 
Fig. 3(6) — Same as Fig. 3(a) with v /V(L) = 10. 
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of the positive from the negative lens. It has the same shape in the hot 
as well as the cold tube and the distance between corresponding points 
of these surfaces in either tube is L, the length of the hot and cold tubes. 
The temperature distribution T/T h in the cold tube is obtained by 
reflecting each point of the temperature distribution of Fig. 3(a) or 
3(b) on the line parallel to the z-axis at T/T h = T h + TJ2T h . The 
neutral surfaces, z/L as a function of x, for various values of flow veloc- 
ity are obtained by rotating the curves of Fig. 4 around the z-axis. At 
high gas velocities (v /V(L) = 10) the neutral surface extends almost 
to the half way point into the hot and cold tubes. 

3.2 Focal Length and Principal Surface 

To calculate effective lenses which describe the ray optics properties 
of the hot and cold tubes it is not permissible to trace rays through each 
tube and compute focal length data from the ray trajectory since each 
tube functions as a combination of positive and negative lenses. It is 
more reasonable to trace rays from one neutral surface to the next 
since the gas between two neutral surfaces acts entirely in one sense 
either as a positive or negative lens. 
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Fig. 4 — Shape of the neutral surfaces for various values of Vo/V(L). 
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We inject a ray at a point z x , x(zi) on the neutral surface with a slope 
[dx/flzfe-zt = and follow it to the point z 2 = Z\ + L, .r(z 2 ). This point 
does not, in general, lie on the next neutral surface since the ray moves 
from its entrance position, x(zi) 9^ .r(z 2 ). However, this point z 2 , x(z 2 ) 
lies sufficiently close to points on the next neutral surface that this ray 
tracing procedure seems justified. 

Our present discussion explains the meaning of the points Z\ and Zi 
(or correspondingly u x and u 2 ) introduced in (13) through (18) and 
shown in Fig. 1. 

The slope and positions of rays entering at z = z x with dx/dz = 
were computed at z = z 2 by numerical integration of (12). The tem- 
perature distribution entering into (12) is given by (19) and (20). 
The values of the slope and the ray position were then used to calculate 
the focal length and principal surface from (14) and (16). The rays 
traveling in the direction opposite to the gas flow were launched at 
2 = z 2 ,x(zi) on the neutral surface with the slope (dx/dz) z=Zi = 
and their slope and position at Z\ = z 2 — L, x(z\) was used to calculate 
P-/L and fJL from (17) and (18). 

It is apparent from (12) and (24) that all of our results depend on a 
parameter 

D = ? -^i T " ~ Tc . (26) 

However, we like to plot our results as functions of W = v /V(L) which 
is contained in a. It is therefore convenient to write 

\> - [J^-A <•['■] (27 



J(L) 
and use 



•© 



L\ . x T h - T r IL 



Cty-in.-U^.-ty (28) 

to characterize the focusing power of the lens.* Fig. 5 shows the focal 
length / divided by the length L of the tubes as a function vq/V (L). The 
solid curves represent the positive, the broken curves the negative lens. 
The focal length of the negative lens is shown as a positive quantity. 
These curves were computed setting x(zi) = 0.1. The positive and 
negative lenses have almost equal focusing power for small values of 
C(L/a). The negative lens has more focusing power for larger values 

* In reference 4 C(L/a) was defined slightly differently. There, T h — T c was 
replaced by Th — '/'■ ('/\ temperature of input gas). 
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Fig. 5 — Normalized focal length f/L of the positive (solid lines) and negative 
(broken lines) lenses of the simple periodic structure as functions of the normal- 
ized flow velocity vo/V(L). 

of this parameter. The temperature gradients are equal but opposite in 
sign for either of these lenses so that one might expect that they should 
have equal focusing power. The discrepancy is explained by considering 
that the rays travel through different parts of the lenses. In the positive 
lens the ray starting at xfa) =0.1 moves closer towards the lens axis 
while the ray in the negative lens, starting at the same point, moves 
away from the axis and toward the wall. 

The minima of the focal length curves are explained by the fact that 
we have no lens action if the gas is stationary, v /V(L) = 0. The lens 
begins to function with increasing gas flow. But, if the gas finally flows 
so fast that the on-axis temperature does not have time to follow, lens 
action ceases again. Interpolation of curves for parameter values other 
than those shown in Fig. 5 is facilitated by noting that the focal length 
is nearly proportional to [C(L/a)]~ l . 

The focal length of an ideal lens does not depend on the input position 
x(zi) of the ray. Plotting f/L as a function of x should result in a straight 
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line parallel to the z-axis. That gas lenses are not ideal lenses is shown 
in Figs. 6(a) through 6(e). These figures contain three different types of 
curves. The solid curves represent the positive lens for rays traveling in 
the same direction as the gas while the dotted curves represent rays 
traveling opposite to the gas flow. The dash-dotted curves give the 
results for the negative lens and rays in the positive gas direction. The 
rays opposite to the gas flow in the negative lens have been ommitted. 
They can be visualized by the fact that the curves in the two directions 
coincide at x = 0. At x = 0.9 the curves for the negative lens join up 
with the dotted lines of the positive lens. The lines for the negative lens 
do not all extend to x = 0.9 because the ray in the negative lens moves 
toward the wall and may hit it before it travels its full length if the lens 
is too strong and if the ray started out sufficiently close to the wall. 

Figs. 6(a) through 6(e) show that there are focal length distortions for 
smaller values of v /V(L). For v /V(L) = 6 and 8 the focal length curves 
are substantially parallel to the .T-axis. (We see, furthermore, that the 
focal length for the two directions of propagation coincide more closely 
for smaller values of C(L/a) and .r. 

The principal surfaces are shown in Figs. 7(a) through 7(e). The 
meaning of the solid, dotted and dash-dotted curves is the same as ex- 
plained above. The dash-dotted lines for the negative lens for the low 
values of C(L/a) coincided very nearly with the solid line for the positive 
lens and was omitted. Also not shown are the corresponding curves for 
the negative lens for rays traveling against the gas flow. The principal 
surfaces are far from being plane. It is also apparent that for most 
values of C(L/a) and x, the two principal surfaces for the two directions 
of the beams don't coincide too closely. This shows not only that the 
lenses comprising the simple periodic structure have considerable dis- 
tortions but also that they are not optically thin under all conditions. 

Fig. 8 shows the dependence of the point x = 0.1 of the principal 
surfaces on the flow velocity v u /V(L). The principal surfaces move to 
z = for vanishing flow velocities and extend far into the tube for large 
values of v /V(L). 

IV. THE EXTENDED PERIODIC STRUCTURE 

4.1 Temperature Distribution 

The extended periodic structure is shown in Fig. 9. It consists of 
alternating positive and negative lenses which are separated by pieces 
of insulating tubes of equal length. 

The temperature distribution, 7' 3 , in the hot or cold tubes is well 
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Fig. G— Normalized focal length f/L as a function of the ray's position x = r/a 
for rays in the positive lens traveling (1.) with the gas stream (solid lines), and 
(2.) against thegas stream (dotted lines) ; and for rays in the negative lens traveling 
with the gas stream (dash-dotted lines). 
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known" ' ' if we can assume that the input gas is at a constant tempera- 
ture. 

T 3 (x,u) = T w + 2(T W - T.) £ /a p x exp (-p H 2 /u) 






(29) 

» > U > W 



T w = either 7\ , temperature of hot tube, or T c , temperature 

of cold tube. 
T i = input temperature to the hot or cold tubes. 
W = v /V(L) with L length of hot or cold tubes. 

The temperature distribution, T t , in the insulating sections is given in 
terms of ^/-functions and their eigenvalues 7 which are defined in the 
appendix. 

T t (xyz) = A 4- 2 A n Un(x) exp -7„ 2 (- - —J . 

W > u > %W. 

The expansion coefficients have to be determined from the condition 

T 3 (x,W) = T 4 (x,W). (31) 

Since the exponential functions appearing in (29) and (30) decrease very 
rapidly with decreasing values of u, it is sufficient to consider only the 
first term in the expansion at the end of each tube. This is justified if 

^ W < 10. (32) 

Condition (31) leads to 

/?o'(l)exp(-/3o 2 /WO 



A = T w - 8(7V - Tt) 



3 fdR \ 



SSa] 



0=00 

and for n j* 

yMd) exp (-tf/W) 



A n = -4(T„- - Ti 



^0(7- - Po > (^--5—) 
where we have used the notation F = dF/dx, 



(33b) 




-**^ 


-^ 






(b) 






i\C> 


H *\ 


cft)= 






\0.75 N 


a75 | 

1 0.05 
V 0.75 




**- = 4 

V(L) 




005 













0.2 












(cj 






v o 

V(L 


= 6 




P^ 


■nsrbj« 


** N 








\ 






eft); 








^ 


\0.3 






\ 
\ 


J.75 


0.05 
\ 0.75 










0.05 




0.2 0.4 0.6 0.8 1.0 0.2 0.4 0.6 0.8 1.0 

1.0 




Fig. 7 — The principal surface for rays in the positive lens traveling (1.) with 
the gas stream (solid lines), and (2.) against the gas stream (dotted lines); and 
for rays in the negative lens traveling with the gas stream (dash -dotted lines). 
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Fig. 8 — The dependence of the point of the principal surface at x = 0.1 on 
the gas velocity. The solid lines represent the positive lens, the dash-dotted lines 
the negative lens for rays traveling with the gas stream. 



The temperature distribution in the hot and cold tubes can be in- 
ferred from curves shown in Ref. 3. The temperature distribution in 
the insulating tubes is shown in Fig. 10(a) for a(L/a) = 0.15 (W = 6.67) 
and in Fig. 10(6) for o{L/a) = 0.05 (W = 20) for various values 
of a[(z — Jj)/a). z — L is the length coordinate counting from the be- 
ginning of the insulating tube. The curves show the temperature as a 
function of radius at different distances from the input to the insulating 



2100 THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 1965 



INSULATION 



HOT 



.__Jf 



COLD 



INSULATION 

i 



INSULATION INSULATION 

I I 

1 COLD r 



Z=0 



z=L 



Z = 2L Z = 3l_ Z = 4L 



Fig. 9 — The hot, cold and insulating tubes comprising the "extended periodic 
structure". 
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Fig. 10(c) — Temperature distribution in the insulating tubes of the extended 
periodic structure following a hot tube. a{(z — L)/a] is the normalized length 
measured from the beginning of the insulating tube. The ratio of input tempera- 
ture to the preceeding hot tube 7\ over its wall temperature T h is Ti/T h = 0.857. 
The normalized length of the tubes is aL/a = 0.15. 

Fig. 10(b) — Same as Fig. 10(a). aL/a = 0.05. The dotted line is the actual tem- 
perature distribution at the end of the hot tube while the line with a [ (z - L)/a] = 
is the temperature distribution at the end of the hot tube which results from 
dropping all but the first term in the series expansion of (29). 
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tube. At z — L — 0, the temperature distribution is identical to that 
at the output end of the hot tube feeding the insulating tube. It is 
apparent that the temperature equalizes rather rapidly. For practical 
purposes we can say that the temperature has reached a constant value 
at <t[(z — L)/a] ^ 0.1. If we consider insulating tubes of length L, 
equal to the length of the hot or cold tubes feeding them, we obtain 
constant output temperatures of the insulating tubes for values ^ 
W < 10. The hot and cold tubes are fed by gas at a constant tempera- 
ture as long as these conditions are met. Therefore, we are justified 
in using (29) which has been derived for the case that the gas at the 
input end of the tube is at a constant temperature. 

In the periodic structure of Fig. 9, the input temperature T to the 
hot and cold tubes are not arbitrary. They adjust themself to satisfy 
the periodicity condition 

(T 3 (x,u)) z=0 = (T 4 (x,u)) z=iL . (34) 

With the help of (34) we can calculate the input temperature T ih 
of the hot or T ic of the cold tube from (29), (30) and (33). 

We obtain from (30) and (33a), for the constant output tempera- 
ture of the insulating tube following the hot tube, (W < 10 is assumed 
so that all exponential terms exp (— y n 2 /W) can be neglected) 

<7v. 0)- = r* - r„ - sm - r„> ft'")ex P (- ft v.n 

and also 



(8Ro\ 
\ ?P A=i 



0o 3 I ~ 



Here we have two equations which allow us to determine the two quanti- 
ties T ih and T ic . 

It is convenient to express them in the form 



T h - T~ 



. /?,/(!) exp j -ftf/W)]- 1 
'BR 



a/3 /,=! 

)3=0o 



(35a) 



2102 THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 1965 



and 



T ie - T c T h - To 



(35b) 



T h - T c T h - T c 
A plot of (35a) as a function of W = v /V(L) is shown in Fig. 11. 

4.2 Focal Length and Principal Surface 

The following graphical representations show the focal length and 
principal surface of one hot and insulating or one cold and insulating 
tube. The extended periodic structure is thus transformed into a system 
of equivalent positive and negative lenses in the same way as in the case 
of the simple periodic structure. Fig. 12 shows the dependence of the 
normalized focal length f/L on the normalized flow velocity v /V(L) for 
a ray entering at r/a = 0.1. The length L is that of the hot tube and not 
the total length of the combination of hot and insulating tubes which 
has the length 2L. The solid curve in Fig. 12 shows the focal length of 
the combination of hot and insulating tubes while the dotted curve 
shows the focal length of the hot tube alone for comparison. It is obvious 
that the insulating tube adds to the focusing power of the gas lens. We 
terminated the curves at v /V(L) = 10 since we wanted to remain in 
the domain of (32) where our simplifying assumptions used to calculate 
the temperature distribution are valid. Fig. 13 shows the corresponding 
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Fj g 11 _ Temperature difference between the hot tube T h and the input tem- 
perature to the hot tube T ih divided by the temperature difference T h - T c 
between hot and cold tube as a function of normalized gas velocity v /V(L). 
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Fig. 12 — Normalized focal length //L of the positive lenses of the extended 
periodic structure (solid lines) and of the hot tubes alone (dotted lines) as func- 
tions of normalized How velocity Vo/V(L). 

curves for the cold lube resulting in a negative lens. A comparison of the 
two figures shows that the negative lens is more powerful than the 
positive lens for corresponding values of C(L/a). 

Figs. 14(a) and (b) show the dependence of the focal length (measured 
from the principal surface) on the input position, x = r/a, of the ray for 
the hot plus insulating tubes. The solid lines represent rays traveling in 
the same direction as the gas while the dotted curves show the focal 
length of rays traveling in Opposite direction to the gas flow. 

The shape of the principal surfaces for the hot plus insulating tube 
are shown in Fig. 15(a) and (b). The distance p of points on the principal 
surface is measured from the gas input end of hot tube. The solid and 
dotted lines again represent rays traveling with and against the gas flow 
respectively. The principal surface is far from being a plane for larger 
values of v Q /V(L). The two principal surfaces do not coincide exactly 
which means that the lens has some optical thickness for larger values 
of C(L/a). 

The corresponding negative lens shows very similar distortions and 
has therefore been omitted. 
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Fig. 13 — Normalized focal length f/L of the negative lenses of the extended 
periodic structure (solid lines) and of the cold tubes alone (dotted lines) as func- 
tions of the normalized flow velocity v v /V(L). 
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Fig. 14 — Dependence of focal length on the ray's input position for the 
positive lenses of the extended periodic structure for rays traveling with the gas 
flow (solid lines) and against the gas flow (dotted lines). 
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Fig. 15 — The principal surface of the positive lenses of the extended periodic 
structure for rays traveling with the gas flow (solid lines) and against the gas 
flow dotted lines. 



V. COMPARISON OF THE TWO PERIODIC STRUCTURES 

In order to compare the two periodic structures let us assume that 
their equivalent lenses are spaced at the same distance D. For the simple 
periodic structure I), the distance between a positive and the next nega- 
tive lens is equal to the length L of the individual tubes. In the extended 
periodic structure D = 2L. It seems fair to compare both structures 
under the condition that the actual gas velocities in either one of them 
are identical. However, this assumption requires some rescaling of the 
data of the extended structure. If we operate the simple structure at 
a certain value of i\,/V(D) the corresponding value for the extended 
structure will be different since v is the same but in the extended struc- 
ture D = 2L. It is apparent that 



Vo 



V{L) 



Vo 



attended 



= 2 



'u 



V{D) ' 



(36) 
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A similar transformation has to be done on C. A certain value of C(D/a) 
of the simple periodic structure corresponds to a value of 



©. 



C - 



= c 



°-) = l -c(° 

2a 4 \a 



(37) 



Since f/L is nearly proportional to C 
convenient to compare the values of 



for small values of C it is 



-m-^m 



This comparison is shown in Fig. 16. For the same values of T h — T c 
and v , the extended structure has the longer focal length because its 
active hot (or cold) tube is only half as long as that of the simple 
structure. The curve of Fig. 16 for the extended structure is not very 
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Fig. 16 — Comparison between the focal lengths of positive lenses of the 
simple and extended periodic structures as functions of the same normalized gas 
velocity v /V(D). 
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accurate for values of v /V(D) > 5 because the value of v /V(D) = 10 
corresponds to Vq/V(L) = 20 for the extended structure. For such a 
large value of the normalized flow velocity, our assumption of a constant 
input temperature to the hot (or cold) tube is incorrect. 

To be able to compare the efficiencies of the two systems we need to 
know the power consumption of one positive lens of either of them 
assuming that it requires no additional power expenditure to keep the 
cold tubes at the temperature T c . 

This power consumption is given by 



P a / 



P = 2irpc p a / {[T(x,u)] z=L - [T(x,u)]^ }xv(x) dx (38) 

v(x) is the gas velocity as a function of x. For viscous, laminar flow 

v(x) = m,(l - x 2 ). (39) 

Using (19) and (24), we obtain for the power consumption per hot tube 
of the simple periodic structure 

f _ 16fl'(l) \ 

irkD(T h -T^ = Vm\ ^ 8 (^ =i (l + exp(/3o7^))| (40) 

and with the help of (29) for the corresponding power consumption 
in the structure composed of extended tubes (assuming that the first 
term in the series of (29) describes the temperature distribution at 
z = I j sufficiently well), we obtain 

2P« t v Tk-T ih ( 1 _ SRp'jl) 



rkD(T k -T e ) V{D)T h -TA tf(aB*\ exp(-j8o7mf. (41) 



\dp/ x 



=i 



The expression (Th — 7' I ; 1 )/(7'/ 1 — T c ) has to be substituted from (35a). 
The quantities represented by (40) and (41) are plotted in Fig. 17 as 
functions of v /V(D). The positive lenses of the extended structure 
consume less power than those of the simple structure for equal amounts 
of gas flowing through them. This is not surprising considering that 
the hot tubes of the extended structure are only half as long as those 
of the simple structure. 

To compare the two structures we require that their lenses have equal 
focal length which we achieve by adjusting the temperature difference be- 
tween the hot and cold tubes. We then plot the ratio of the resulting 
power consumptions. This plot is shown in Fig. 18. The extended struc- 
ture requires less power than the simple one and this ratio improves as 



2108 THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 1965 















SIMPL 
STRUCTU 


RE/^ 








/ 


EXTENDED 
STRUCTURE,, 

























Fig. 17 — Comparison between the power consumptions of the positive lenses 
of the simple and extended periodic structures as fund ions of Ihe same normal- 
ized gas velocities v»/V(D). 
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Fig. 18 — The ratio of the power consumption of the positive lenses of the 
extended periodic structure over the power consumption of the positive lenses 
of the simple periodic structure for equal focal length as a function of normalized 
gas velocity vo/V(D), 
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the flow velocity i ncre a ses. It should be remembered, however, that the 
focal length of the extended structure (Fig. 16) is inaccurately repre- 
sented for values v /V(D) > 5. This inaccuracy carries over to Fig. 18. 
Nevertheless, it is clear that the efficiency of the extended structure is 
more than 20 per cent higher than that of the simple structure for 
v /V(D) > G. The additional focusing which the insulating tube sections 
provide pays off to some extent to make the extended structure more 
efficient. The improvement is not as large, however, as one might have 
hoped and it may be considerably poorer if the insulating tubes have the 
finite heat conductivity of a real material. 
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APPENDIX 

The functions to be discussed in the appendix are solutions of the 
differential equation 

^ + 1^ + ^(1 - /)F = 0. (42) 

ax- x dx 

The solutions of (42) are related to Whittaker's functions W tJt by 

F = -TF a . /4 ).o(fc.r 2 ). 
x 

The differential equation (42 ) stems from an approximate formulation 
of a heat transfer problem. 7 Assume that a gas at a different tempera- 
ture is blown into a tube with a circular cross section. The gas is sup- 
posed to flow as a viscous fluid in laminar flow with a velocity profile 

v(r) = v (l -^J (43) 

/• = distance from tube axis 
a = tube radius. 

The stationary state of the temperature distribution T is obtained from 

a V 2 T = WT (44) 

with a being a constant which contains the heat conductivity, density 



2110 THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 1965 

and specific heat at constant pressure — all of which are assumed to 
be temperature independent which, strictly speaking, is not true. 

The velocity has only a longitudinal component v z given by (43). 
In polar coordinates r, <p, z taking d/d<p = (44) can be written 

jd 2 T , ldT , d 2 r\ A Aar ,._* 

a \W + 7JF + l*I- H \ 1 -*)*' (45) 

It is often permissible to neglect d 2 T/dz 2 compared to the term on the 
right hand side of (45). Taking 

r = ax (46) 

and 

T(r,z) = Fix) exp (-Xz) (47) 

we get from (45) 



■1 



Finally, taking 



** + lf + *!* (!-*> = 0. (48) 

ax- x ax a 



X = k 2 4 (49) 

a 2 v 



we recognize that (48) is identical to (42). 

We are interested in two types of heat transfer problems : 

(1.) The tube through which the gas flows may be kept at a constant 

temperature T w . We then have the boundary condition T{a,z) = T w . 

It is more convenient to introduce a new variable 

d(r,z) = T w - T(r,z). (50) 

We can replace T by 8 without changing any of equations (44) through 
(47 ) . However, the boundary condition now becomes simply 

6(a,z) = 0. (51) 

The functions satisfying this boundary condition are designated by 

F(x) = B(x) (52) 

and (51) becomes 

72(1) = 0. (53) 

The R functions have been studied in some detail. ' 

(2.) The second type of problem involves a tube which is a perfect 
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heat insulator. That means that no heat flows into or out of the walls. 
This assumption requires that 



m 



, =0. (54) 

Or /r=a 

A second class of functions is obtained by setting F(x) = U(x). The 
[/-functions satisfy the same differential equation but are defined by 
the boundary condition 



WA=i 



^ =0. (55) 



For convenience, both functions are normalized so that 

22(0) = f/(0) = 1. (56) 

The [/-functions have not been studied to the knowledge of the author. 

For series expansions of arbitrary heat distributions in terms of either 
the R or the U functions we need their orthogonality relations and any 
numerical evaluation of heat transfer problems requires the knowledge 
of the eigenvalues and numerical values of these functions. 

The eigenvalues belonging to the R functions R n (x) will be designated 
as 

h = A. (57) 

and those of U n (x) will be designated by 

fc n = 7» • (58) 

Orthogonality Relations 

Let F n with eigenvalue k„ designate either an R„ or a U„ function. 
We proceed to show that 



'0 

We have 



and 



f x(l - x 2 )F n (x)F m (x) dx = for n 5* m. (59) 

Jo 

F n " + - Fj + K\l - x)F n = (60) 

X 

F m " + - F-J + fc m 2 (l - x 2 )F m = 0. (61) 
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We multiply (60) by xF m and (61) by xF„ , subtract and integrate: 

(k n 2 - kj) f x(l - x 2 )F n (x)F m (x) dx 
Jo 

- - £ ■ {{"- ( f »" + i fJ )) - f ° ( f "" + i F -')} dx - 

We perform partial integrations and obtain for the right hand side of 
this equation 

-[xF m F n ' - xF n F m % + f \F m F n ' - F n F m ' + x(F m 'F n ' - F a 'F m ') 

Jo 

- (F m Fn - F n F m ')\ dx = (F m F n ' - F n F m ') x=l = 0. 

The last part of the equation follows from the boundary condition (53 ) 
or (55), depending whether F stands for an R or a U function. This 
calculation proves (59). 

Next, we calculate the value of (59) in the case n — m. 

[ ,(1 - *>)[FM? dx = fcn fig I F k f \ =i 

X 2fc„ h=\ n - 

For F» = R n we get 

f x{\ - x 2 )[R n (x)f dx = -i- [|f fi/L (62) 

Jo 2/3„ \_d(3 J0=ft, 

and for /''„ = U n 

f x(l - x 2 )[U n (x)f dx = -±- \u n A d/L • (63) 

Finally, we need to determine the value of the integral over the prod- 
ucts R„U m : 

Rn" + - Rn + /3n 2 U - Z 2 )# n = (64) 

X 

U m " + i Uj + 7n 2 (l - x 2 )J7 m = 0, (65) 
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(/3„ 2 - Tm 2 ) f 3(1 - x)R»(x)U m (x)dx 
Jo 

= [ {x(R n U m " - U m Rn") + R n U m ' - U m Rn]dx 

J I) 

= [RnUj - U m R n '] x -l 



or using (53) and (55), 

f x(l - af )«„(*) tf.UKr = , 1 , [0»ti . (66) 

Calcvlation of the R and U Functions and Their Eigenvalues 
We make the series expansion 

F(x) = ZC 2 ,T 2 ". (67) 

f=0 

For the problem of interest to us F(x) has to be an even function of x, 
for that reason only even powers of x appear in (67). The normalization 

F(0) = 1 

requires that (68) 

Co = 1. 

The substitution of (67) into (42), using (68), leads to 



and 

C 2e = -^- {C 2B _4 - C 2E _ 2 } for v ^ 2. (69) 

(2w) 2 

The parameter k has to be chosen so that either F(l) = Oor^'(l) = 
results, depending whether k and F shall represent /3 and R or 7 and U 
respectively. 

The fact that k 2 enters all coefficients C- iv makes the determination 
of /3„ and y n very tedious. 

A further difficulty results from the fact that the coefficients C%, 
grow to very large values particularly for the larger values of /3„ and 
y„ before they decrease again. The series (67) does not converge readily 
for values of x close to 1 . In fact, it proved impossible to compute more 
than the first eight R and U functions from (67) on the IBM 7094 com- 
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puter even using double precision since the absolute value of R and U 
remains between zero and one but the coefficients C 2 „ grow to values 
above 10 20 . The series (67) can be used to compute R n and U n for x 
in the range ^ x ^ 0.5 since the powers of x decrease rapidly enough 
to keep the value of the product Cz v x 2v within manageable proportions. 
However, in order to cover the whole range ^ x ^ 1 it proved 
necessary to use the following series expansion: 

00 

F(x) = Z A/ with y = 1 - x (70) 

i>=0 

to calculate R and U in the range 0.5 ^ x ^ 1. 
Equation (42) expressed in terms of y reads 

d - v) ^ - f + k2 ^y - *v + y 3 ^ F = °- (71 > 

dy 2 dy 

The coefficients A and A have to be properly chosen to satisfy the 
boundary conditions at x = 1. For F = R we require R(\) = so 
that 

Do = 

results. For F = U we require £/' (1 ) = so that 

A = 

results. 

The substitution of (70) into (71 ) yields 

D 2 = iA , A = £A - ^- 2 A , A = (1 - £fc 2 )A 

1 (72) 

A = 7 7t {(» - D 2 -D,-i - fc 2 (2A- 3 - 3A-4 + A-s)}. 

y(y — 1) 

The eigenvalue k = p or k = y and the coefficient A or A must be 
chosen so that F as well as F' are continuous at x = 0.5 where both series 
expansions should coincide. 

By breaking the range of x into two parts and using different series 
expansions to cover both parts of the range it was possible to compute 
the function and their eigenvalues. Table I shows the eigenvalues 
/3„ and y n as well as the values of dR„/d@, — A = Rn- (d/dy)U n ' and 
Do = U„ all taken at x = 1. These values are needed to evaluate the 
integrals (62), (63) and (66). 

The values of dR n /d(3 and dU n '/dy were obtained from differentiation 
of the series (67) and evaluating it at x = 1. The terms of the differ- 
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entiated series grow very large so that only the first eight values of 
dR/d(3 and the first six values of dU'/dy could be obtained. The remain- 
ing values of dR/d(3 were calculated from the approximation 8 



d_R, =(_!)» 



d/3 /*=! 



6 2 ' 3 r (§) /v /3 



(73) 



which is in good agreement with the values obtained by machine calcu- 
lation for larger values of n. An approximation of dll/dy can be obtained 
by using approximations similar to the ones used for the R functions 
in Ref. 8. One gets 



7=7 



™ =-(-!)" 



7T7 



6sr (i) 



(74) 



However, this approximation is not very good for n ^ 15 so that we 
used the equation 

" A u 



(f'L --<-»■* 



7T7 



6 

m=i 7" 



(75) 



The coefficients A^ were determined from the first six values of dll'/dy 
which were obtained from a machine calculation. Their values are given 
at the bottom of Table I. 



Table I 







Rn (1) 


(9RW\ 

\ 3(3 / p = p n 


In 


tf«(U 


\ Sy Jy = y n 









2.7043G 


-1.01430 


-0.50090 





1 




1 


6.67903 


1.34924 


0.37146 


5.06750 


-0.492517 


0.97816 


2 


10.67338 


-1.57232 


-0.31826 


9.15760 


0.395509 


-1.24720 


3 


14.6711 


1.74600 


0.28648 


13.1972 


-0.345874 


1.43522 


4 


18.6699 


-1.89090 


-0.26449 


17.2202 


0.314047 


-1.58486 


5 


22.6691 


2.01647 


0.24799 


21.2355 


-0.291252 


1.71127 


6 


26.6686 


-2.12814 


-0.23491 


25.2465 


0.273806 


-1.82164 


7 


30.6682 


2.22038 


0.22485 


29.2549 


-0.259853 


1.92042 


8 


34.6679 


-2.32214 


-0.21548 


33.2615 


0.248332 


-2.01037 


9 


38.6676 


2.40274 


0.20779 


37.2669 


-0.238591 


2.09330 


10 


42.6667 


-2.48992 


-0.20108 


41.2714 


0.230199 


-2.17045 


11 


46.6667 


2.56223 


0.19516 


45.2752 


-0.222863 


2.24275 


12 


50.6667 


-2.64962 


-0.18988 


49.2785 


0.216371 


-2.31088 


13 


54.6667 


2.70216 


0.18513 


53.2813 


-0.210569 


2.37539 


14 


58.6667 


-2.76421 


-0.18083 


57.2837 


0.205216 


-2.43671 


- 


4. = 


-4.881355 


.4 4 = 2.838701 10 4 






A 2 = 


1.536461 10 2 


A B = -1.420240 10 6 






A 3 = 


-2.838383 10 3 


.4 6 = 2.728875 10 s 
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All approximate formula for /3 is s 

fti -4n+ |. (76) 

The /? values of Table I f or n ^ 11 have been computed from (76). 
The lowest order U function is a constant 

U = 1 with to = 0. (77) 

The 7 values should converge to 

7n = 4n + | . (78) 

This expression can be derived by methods analogous to those used in 
Ref. 8. For unknown reasons this approximation is much poorer than 
that for |8 n . However, it appears from the values of Table I that y n 
will converge to (78 ) for very high values of n. 
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